Monomer-Dimer Mixture on a Honeycomb Lattice 
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We study a monomer- dimer mixture defined on a honeycomb lattice as a toy model for the spin 
ice system in a magnetic field. In a low-doping region of monomers, the effective description of this 
system is given by the dual sine-Gordon model. In intermediate- and strong-doping regions, the Potts 
lattice gas theory can be employed. Synthesizing these results, we construct a renormalization-group 
flow diagram, which includes the stable and unstable fixed points corresponding to Mz, and Me in 
the minimal models of the conformal field theory. We perform numerical transfer-matrix calculations 
to determine a global phase diagram and also to proffer evidence to check our prediction. 

PACS numbers: 75.40.Cx, 05.50.-|-q, 05.70.Jk 



An introductory study of dimer degrees of freedom for 
condensed matter physics was carried out long time ago. 
In a paper by Fowler and Rushbrooke in 1937, the dimer 
represented a diatomic molecule adsorbed on a crystal 
surface [ij. Later, the statistical mechanics of dimers was 
studied by Kasteleyn and Temperley and Fisher Q. 
These pioneering works have provided exact solutions of 
two-dimensional models under certain conditions. Since 
then, dimer models have gathered increasing attention. 

The correspondence between dimer models and real 
substances seems to fall into two categories. First, like 
for the case of diatomic molecules, the correspondence 
originates in the shape of the elements. For example, re- 
cently a network system made of rodlike molecules on a 
substrate was measured by using scanning tunneling mi- 
croscopy It was argued that, because of the shape of 
the molecules, a honeycomb-lattice dimer model can be 
employed to reproduce properties of the network . Sec- 
ond, the correspondence can originate in interactions be- 
tween elements. Spin configurations in the ground state 
of a triangular-lattice antiferromagnetic Ising model can 
be related to those of dimers on a honeycomb lattice 0] . 
In this case, dimer degrees of freedom emerge on a dif- 
ferent lattice and represent unsatisfied bonds. 

In this work, we investigate a dimer-based model which 
can be related, in the latter sense, to the spin ice in a 
magnetic field. To describe their relevance, we begin by 
summarizing related researches. The rare-earth titanates 
such as i?2Ti207 {R = Ho, Dy) are known as the Ising py- 
rochlore magnets (IPM) [7|, where each R^^ behaves as 
an Ising spin along a local axis pointing in the center of a 
tetrahedron. In the low temperature, due to the magnetic 
couplings 0, six states, satisfying the so-called 2-in-2-out 
condition (also known as the ice rule), are permitted for 
each tetrahedron. When a magnetic field is applied along 
its [111] direction, a plateau is observed at | of the sat- 
uration magnetization [l^ . This incompressible state 
is called the kagome ice. For a kagome layer sandwiched 
between two triangular layers, only three states are per- 
mitted for each tetrahedron pointing in the [111] ([111]) 
direction since spins on the triangular layers are polar- 



ized. Then, the mapping from the spin configuration on 
the kagome layer to the dimer configuration on the hon- 
eycomb lattice is possible, which enables us to exactly 
enumerate the residual entropy [l]| . 

While the magnetic field causes a dimensional reduc- 
tion, it also induces low-energy excitations that lead to 
breakings of the ice rule: For each tetrahedron pointing 
in the [111] ([111]) direction, a state that minimizes the 
Zeeman energy takes the 3-in-l-out (l-in-3-out) configu- 
ration. This state can be viewed as a magnetic monopole 
(antimonopole) in the solid [l2-15|. Specific heat mea- 
surement revealed that the temperature dependence of 
the monopole density obeys the Arrhenius law with an 
energy scale controlled by the magnetic field 14|. So, 



monopole excitation effects become significant at least 
around the phase boundary between the kagome ice and 
a saturated state [l2,|T5|. We note that since monopole- 
antimonopolc pair creations and annihilations are only 
permitted, it follows naturally to introduce monomers in 
the dimer model to describe them. 

Motivated by the research advances on IPM, we con- 
sider a monomer-dimer mixture (MDM) on a honeycomb 
lattice Afj . When we write the dimer occupation number 
on the kagome lattice Ak (the medial lattice of Ah) as 
7T,; = 0, 1, then a reduced Hamiltonian is given by 



p=a,b,c {l,m)£A^ 



ni 



(1) 



^•^^b,c^ three sublattices of Ak). The first term repre- 
sents the interaction between two neighboring dimers and 
the second term signifies the chemical potential that con- 
trols the density of dimers (and thus monomers). In ad- 
dition, we impose the hard-core constraint, Vj S Ah, 
~ 1' where {l{j)} denotes three sites 
around j. For the compounds under consideration, the 
long-range dipole interaction is expected to be large Q. 
While the so-called projective equivalence explains their 
spin-ice behaviors [l6|, the long-range nature may be- 
come important with breaking the ice rule. Here, we 
have employed the oversimplified model with the short- 
range interaction because we focus on universal proper- 
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ties stemming from a competition between the interac- 
tion and the monomer doping effects, which could shed 
some light on the understanding of, for instance, the 
properties of monopoles in the kagome ice. 

The MDM system on the square lattice was discussed 
by several authors 
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[18|, where an effective theory 
for long-distance behaviors plays an important role. We 
take a similar approach; however our analysis predicts 
the emergence of criticalities not observed in the square 
lattice case. Let us start with the dimer covering case 
without the interaction (u = 0, /i = oo). Its critical be- 
haviors are described by a two-dimensional sine-Gordon 
Lagrangian density = Cq 



Cip with 



cos 



3V2( 



(2) 



We have denoted the coarse-grained height field as (p{x.), 
which satisfies the periodicity condition ^/2(p = V2ip + 
2ttN with iV e Z (a, a short distance cutoff) [itI. Ii9|. 
Since the Gaussian coupling K equals i, the nonlinear 
term that represents the discreteness in the original 
height field is irrelevant (y^, a negative constant). 

Next, we consider a modification of Cq that is brought 
about by the nearest-neighbor interaction. Jacobsen and 
Alet discussed the same effect in a somewhat different 
context [1] and concluded that, like the square-lattice 
case, it is devoted to a renormalization of the Gaussian 
coupling [13, [2^ • Thus, we can write K{u) ~ i+ciw (ci, 
a negative constant). On the other hand, a monomer on, 
say, the A (B) sublattice of Ah corresponds to a defect 
with a positive (negative) charge. In the sine-Gordon 
language, using the disorder field 1? dual to (p, we can 
express monomers as e*'^''. Consequently, at least in 
the low-doping region, a dual sine-Gordon model C = 
^g{u) + C^p + describes the MDM system, where 

cos V2^ : (3) 



1? 
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and a fugacity of monomers cx e"''/^. To proceed 
further, we clarify a renormalization-group (RG) flow in 
the attractive region u < 0. Since the dimensions of 
cosSv^"^ and cosv^'? are 9/2K and K/2, respectively, 
the Gaussian fixed line is realized for i < X < | in 
the dimer covering case, whereas it is unstable against 
the doping. We sketch the flow in Fig. [1] (see the red 
region), where (w, z) = (e~", e"''/^). The open circle 
on the w axis denotes the end point K ~ |, where 
Cip becomes marginal and brings about the Berezinskii- 
Kosterlitz-Thouless (BKT) transition to the 3-fold de- 
generate columnar ordered state, as shown by the flow 
to the filled triangle. Meanwhile the relevant £^ leads to 
the disordered phase (see the flow to the filled circle) . dp 
and £^ are mutually non-local, so intrinsically they are 
in competition. If both of them are relevant, it becomes 
obvious as a flow to an infrared (IR) fixed point. Since 
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FIG. 1: (Color online) A rough sketch of the RG flow diagram. 
The open circle denotes the end point of the Gaussian fixed 
line on the w axis; the marks for fixed points are given with 
explanations. The sine-Gordon theory (the Potts lattice gas 
theory) is varied in the red (blue) region. 



this fixed point describes the transition to the 3-fold de- 
generate ordered state, it may be of the 3-state Potts 
universality class. Consequently, we can expect a mass- 
less RG flow from the ultraviolet (UV) fixed point with 
the central charge c = 1 (the open circle) to the IR fixed 
point with c = | (the open triangle) [21j. Indeed, C is 
an effective theory for the 3-state Potts model, and, un- 
der the self-dual condition, this flow has been explicitly 
obtained and understood as a renormalization of the Z3 
neutral X operator, with dimension xx = "y, to the IR 
fixed point [2l[ . Although C is not self-dual around the 
open circle, we suppose that the same RG fiow should be 
observed at least around the open triangle. We provide 
numerical evidence for this assumption below. 

For the exploration of intermediate- and strong-doping 
regions, let us focus on the role of the X operator. Writ- 
ing the minimal model series of the conformal field the- 



ory with c = 1 — 6/p{p + 1) as Aip and a primary field 



at the position (r, s) of the conformal grid as 



a(p) 



then 



X 



^(5) 



i.e., the leading irrelevant operator on 
Zamolodchikov [22| and Ludwig and Cardy Q discussed 
the deformation of A4p by the least relevant operator (jf^^l 
and concluded that there exists a RG fiow connecting 
minimal models as 



UV : M 



IR : 7Wp_i, 



'^3,1 ' 



(4) 



where a UV-IR operator transmutation is also given. By 
taking these into account, it is plausible that another UV 
fixed point exists and that a RG fiow connects Aie (the 
open square) and Ai^ (the open triangle), as given in Fig. 
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[TJ In the opposite direction, the renormahzation always 
amphfies a deviation from Aie, so the coupling of 
flows to a strong-coupling fixed point (the filled square). 
It may be a discontinuity one, and thus the transition is 
of first order in the strong-doping region. 

The emergence of the new critical fixed point Mq with 
c = I is key to understanding the MDM system. We take 
a look at the following UV-IR operator correspondence 
and characterize it 1231 . 
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UV : 0; 



4% 



IR 



^(5) 
"3,1' 



^(5) 
''2,1' 



(5) 



The latter two of the three IR operators correspond to the 
energy density (e) and the Z3 spin (a) of the 3-state Potts 
model, respectively. Thus, (= as well as ^^^3 

is Z3 neutral, and represents a thermal operator in the 
Lagrangian. Since, using the Kac formula, dimensions 
of UV operators are given as x^^l = ^, x'fl = |, and 

Xg^g = ^ , ^^^3 and are both relevant and provide the 
flow around the open square in Fig. [TJ In the literature, 
one can find a model that exhibits the same RG flow 
given in the blue region of Fig. [1] Nienhuis discussed the 
g-statc Potts lattice gas model on a square lattice [23 |. 
where the leading and next-leading thermal exponents 
were obtained on both the critical and tricritical points. 
Indeed, for (7 = 3, those on the latter are yi = | and 

y2 = which agree with the dimensions x^^g and x^^\. 
Therefore, the MDM system in the intermediate- and 
strong-doping regions can be viewed as the 3-state Potts 
lattice gas, where monomers and dimers play roles of the 
vacancies and the Z3 spins, respectively. 

Obviously, the 3-fold axis symmetry in Ajj plays a de- 
cisive role in determining Cip and mapping to the Potts 
lattice gas with g = 3. In contrast, for the MDM on 
the square lattice [H, [11], due to the 4-fold axis sym- 
metry, the sine-Gordon theory with a potential cos4-\/2v' 
and a 4-state Potts lattice gas theory are relevant in the 
weak- and strong-doping regions, respectively. So, in- 
stead of crossovers of criticalities, we observe a fixed line 
that terminates by the first-order phase transition line. 

Now, we explain numerical calculations and results 
to check our prediction. We have focused on a struc- 
ture of the phase diagram and the universality classes of 
phase transitions. We first summarize our results on the 
phase diagram in Fig. [2j We quote the BKT point of Q , 
Wbkt — 1/0.635, which is indicated by the open circle on 
the w axis. The curve for open triangles (filled squares) 
denotes the second-order (first-order) phase transition 
boundary Z2{w) [zi{w)] between the ordered and disor- 
dered phases. The open square represents the tricriti- 
cal point. For their enumerations, we have performed 
numerical transfer-matrix calculations of Eq. ([T]) on Ah 
with L X 00 cylinder geometry [L, the circumference of 
the cylinder in units of the length between two parallel 
bonds, is given in multiples of 3 because the columnar 
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FIG. 2: (Color online) The global phase diagram. Open tri- 
angles (filled squares) represent the second-order (first-order) 
phase transition points; the open square denotes the tricritical 
point. The dotted line gives the asymptotic form z — w. 



state is expected]. We note that, with the monomer dop- 
ing, the row-to-row transfer matrix T(L) becomes less 
sparse. Thus, the accessible system size is strongly lim- 
ited to a small number, e.g., L < 12 in our calculations. 

To determine the second-order phase transition points, 
we have performed phenomenological RG (PRG) calcu- 
lations. Writing the eigenvalues of T(i) as Xi{L) and 
their logarithms as Ei{L) = —\n\Xi{L)\ {i specifies an 
excitation), then the conformal invariance provides the 
direct expressions of c and Xi in the critical systems as 
Eg{L) ~ Lf ~ 7rc/6iC and A£;j(i) ~ 2nx^/LC 
Eg{L), AE,{L) [= E,{L) - E^{L% C (= a/3/2), and / 
correspond to the ground-state energy, an excitation gap, 
the geometric factor, and the free-energy density, respec- 
tively. Eg is found in the zero momentum sector (fc = 0). 
Also, Ei is in a sector specified by its symmetry property. 
In the PRG calculation, we have employed the lowest ex- 
citation, corresponding to cr, with k = 27r/3. We then 
numerically solve the condition LAEa-{L) = L' AEa-{L') 
with respect to z for given values of w. We have ex- 
trapolated finite-size estimates of L = 6, 9, and 12 to 
the thermodynamic limit based on an assumption of the 
leading 0{L~^) correction [27[. Then, we find that the 
curve Z2{'w) continues smoothly to the open circle i«bkt, 
which exhibits a consistency with the previous result (Hj . 

To check the universality class, we evaluate c and Xa 
along Z2{w). Plotted in Figs. [3l[a) and[3l^b) are the ex- 
trapolated results [i^ . Although we cannot extract reli- 
able data around i«bkt due to the smallness of the sys- 
tem size, we can find evidence to support the 3-state 
Potts universality c = | and x„ — down to u ~ —1.5. 
Also, Fig. [3lc) provides a finite-size scaling plot of the 
gap AEa(L) = - Z2{w)]L^} aX u = -1. Since 
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FIG. 3: (Color online) The estimations (open circles) of (a) c 
and (b) Xa- along Z2{w). Solid and dotted lines denote theo- 
retical values for Ms and A4e; the crosses also give those at 
the BKT point. The vertical arrow indicates the tricritical 
point, (c) The finite-size-scaling plot of AEa at u — —1. 



Dy2Ti2 07 revealed a first-order phase transition between 
the kagome ice and the saturated state at 7J ~ 0.9 T for 
T < 0.36 K In contrast, in our MDM system, the 

critical phase is not stabilized in the doped region. This 
implies that the long-range dipole interaction plays a cru- 
cial role for a stabilization of the kagome ice observed in 
real materials. However, this scenario is currently at the 
level of speculation, and its confirmation is left as a future 
work. 

The author thanks K. Goto, H. Kadowaki, G. Tatara, 
M. Fujimoto, and K. Nomura for stimulating discussions. 
Main computations were performed by using the facilities 
at the Cyberscience Center in Tohoku University. 



z — Z2 linearly couples to e, we expect the plot with 
ly ~ 1/(2 — .Te) = I to yield a collapse of the finite- 
size data onto a single curve. The result exhibits a good 
scaling property and thus supports our prediction. On 
one hand; for u < —1.5, we find deviations of data from 
the values, which implies that the system is approach- 
ing the tricritical point. We search it along Z2{w), using 

„(6) 



the criterion 3 = ^ [see Fig. [S^b)] [28| and estimate 
{wt, zt) ~ (5.38, 5.03), which is given by the open square 
in Fig. [5] (the monomer density is estimated as around 
25%). To check its criticality, we estimate the central 
charge and scaling dimensions of the thermal operators. 
The results are c ~ 0.859, x[^l ~ 0.293, and xj*^] ~ 1.439. 
We thus find reasonable agreement (within a few percent) 
between the numerical data and the theoretical predic- 
tions which strongly supports our RG argument. 

In the strong-doping region, we have found some data 
to imply a first-order phase transition, e.g., an abrupt 
change in the monomer density and a double-peak struc- 
ture in the monomer-density distribution function. These 
could be used as means to determine the phase transition 
boundary zi{w). However, we provide here numerical so- 
lutions of the PRG equation, represented by filled squares 
with the curve in Fig. [2] because they are known to give 
accurate estimations |29|. For large z and w, we can 
analytically estimate the phase boundary via an energy 
comparison between the complete columnar state and the 
dimer vacuum. For the honeycomb lattice case, z = w, 
which is given by the dotted line in Fig. [2] We observe 
that filled squares asymptotically converge to the line. 

In conclusion, motivated by recent experiments on 
IPM, we have investigated the MDM system on the hon- 
eycomb lattice. We have provided the global RG flow 
diagram, which determines universality classes of phase 
transitions. Also, we have performed numerical calcu- 
lations and provided evidence that supports our predic- 
tions. Regarding the relationship to experiment, there 
exist some issues. For instance, the phase diagram of 
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